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ABSTRACT 

The Apache Point Observatory Galactic Evolution Experiment has measured the stellar velocities 
of red giant stars in the inner Milky Way. We confirm that the line of sight velocity distributions 
(LOSVDs) in the mid-plane exhibit a second peak at high velocities, whereas those at |6| = 2° do 
not. We use a high resolution simulation of a barred galaxy, which crucially includes gas and star 
formation, to guide our interpretation of the LOSVDs. We show that the data are fully consistent 
with the presence of a thin, rapidly rotating, nuclear disk extending to ~ I kpc. This nuclear disk is 
orientated perpendicular to the bar and is likely to be composed of stars on x2 orbits. The gas in the 
simulation is able to fall onto such orbits, leading to stars populating an orthogonal disk. 

Subject headings: Galaxy: bulge — Galaxy: evolution — Galaxy: formation — Galaxy: kinematics 
and dynamics — Galaxy: stellar content 


1. INTRODUCTION 

Detections of high Galactic standard-of-rest velocity 
(Lgsr) peaks in the Apache Point Observatory Galactic 
Evolution Exper iment (APOGEE) co mmissioning data 
were reported bv iNidever et al.l (I20I2I1 across all fields at 
4° < I < 14° and —2° < b < 2 . iNidever et ahl (|2012[ 1 
interpreted the high-VcSR peaks as being due to stars in 
the Galactic bar. However, the pea ks are not sta tisti- 
cally significant in a number of fields ()Li et al.ll20I4fl and 
no high-VcSR peaks were found at negative longitudes in 
the Bulge Radial Ve locity Assay (BRAVA), at b —4° 
(iKunder et al.l[2Mll . Additionally, no h igh-UnsR peak s 
can be found in p ure iV-body models ()Li et al.l 120141 1. 
iMollov et al.l (|2015[ 1 demonstrated that resonant (2:1 and 
higher order) orbits, viewed on their own, were a ble to 
generate high-VcSR peaks. lAiimer fc Schonrichl (l2015f l 
proposed that such resonant orbits are populated by 
young stars recently trapped by the bar; they argued 
that the APOGEE selection function is biased toward 
such young stars. 

Bars have been implicated in building large gas reser¬ 
voirs at the centers of galaxies, fuelling high star forma¬ 
tion rates there. As in other barred galaxies, the Milky 
Way (MW)’s bar funnels gas inwards (IBinnev et al.l[T991l : 
iWeiner fc SellwoodI 119991 : iFuxI 1199911 . This gas gives 
rise to structures such as the Gentral Molecular Zone 
(GMZ), spanning —1° < I < 1.5°. The GMZ contains 
5 —IO xIO'^Mq of molecular gas (|Ballv et al.lll987HGustenl 
1198911 ■ driving a star formati on rate of ^ 0.14 MQyr“^ 
( Wardle fc Yusef-Zadehll2008ll . A molec ular gas disk ex¬ 
tends across |^| < 6° and |6| < 1.6° (iBovce &: CohenI 
11994 [Pa me fc T haddeusl[l994ll . iLiszt fc BurtonI (Il980f l 
and iFerriere et al.l (j2007f l interpreted the observed molec¬ 
ular, atomic and ionized gas outside the GMZ to Galactic 
longitude \l\ ^ 10° as a (tilted) disk with semi-major axis 
of radius ~ 1.4 kpc with a hole at its center. In external 
galaxies, star formation in nuclear rings builds nuclear 
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disks (jKormendv fc KennicuttI 1200411 . In this Letter we 
demonstrate that the high-VcSR peaks in the line of sight 
velocity distributions (LOSVDs) are consistent with the 
presence of a nuclear disk in the MW. 

2. SIMULATION 

Here we use a high resolution simulation, with gas 
and star formation, which develops a bar, dri ving gas to 
the c enter and forming a stellar nuclear disk (iGole et al.l 
l2014f) . to derive the kinematic signatures of such a disk. 
We use these to g uide our interpreta tion of the APOGEE 
Data Release 12 (jAlam et al.ll2015ll stellar velocity data 
for the inner MW. While the s i mulat ion was not designed 
to match the MW, IGole et Ml (j2ni4[ l showed that the nu¬ 
clear disk that it forms is qualitatively similar to those 
in external galaxies. 

The simulation was evolved with the N- 
body-l-smoothed part icle hydrodynamics code gasoline 
(|Wadslev et al.l 120041 1. The galaxy forms out of gas 
cooling off a hot corona in pressure equilibrium within 
a dark matter halo of virial mass M 200 = 9 x 10^^ Mq. 
Both the dark matter halo and the initial gas corona are 
represented by 5 x 10® particles. As the gas cools and 
reaches high density, star formation is triggered. Star 
particles then provide feedback via winds from massive 
stars, and types la and H supernovae (iStinson et al.l 
[ 2 OOI). Gas particles all have initial mass of 2.7 x IO'^Mq 
and star particles are spawned from gas with 35% of 
this mass. This high mass resolution allows us to use a 
high star formation th reshold of 100 cm“® for the gas 
(iGovernato et al.l 12010(1 . By the end of the simulation 
the galaxy has a stellar mass of 6.5 x 10^® Mq in 
^ 1.1 X 10' particles. This large number of star particles 
provides a very fine sampling of the mass distribution 
at the center of the mo del. Further de tails of the 
simulation are provided in iCole et all (j20l4l 

The bar forms at around 4 Gyr. After 6 Gyr a promi¬ 
nent nuclear disk starts to form which, by 10 Gyr, has 
a semi-major axis of 1.5 kpc. The nuclear disk is per¬ 
pendicular to the bar and its stellar streaming is perpen¬ 
dicular to the bar’s. At 10 Gyr the nuclear disk in the 
simulation is quite massive and is thus unlikely to match 
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any nuclear disk in the MW. Therefore here we consider 
the model at two earlier times: at = 6 Gyr, before the 
nuclear disk forms, and at ^2 = 7.5 Gyr when a strong 
nuclear disk is established. Aside from the nuclear disk 
becoming more massive and the bar growing longer, the 
model at 10 Gyr is not qualitatively different from at ^ 2 - 

2.1. Scaling to the MW and Viewing Perspective 

In order to compare to the MW, we rescale the model in 
both size and velocity. Size rescaling is accomplished by 
matching the size of the bar to that of the MW. Between 
ti and ^2 the average size of the bar in the simulation, 
as measured from the radius at which the phase of the 
m = 2 Fo urier moment deviat es from a constant by more 
than 10° (lAguerri et al.ll200^ . is 2.1 kpc. Assu ming that 
the M W’s bar has a semi-major axis of 3.5 kpc (|Gerhardl 
[ 2 ?iol . we scale all coordinates by a f actor of 1.67 . (Sca l- 
ing to the more up-to-date bar size of lWegg et al.l (j2015[ l. 
5 kpc, leads to a nuclear disk which is much too large; 
because we seek a closer nuclear disk size match, we scale 
to the older bar size, but this is not to imply that the 
real MW bar semi-major axis is closer to 3.5 kpc than 
5 kpc.) The velocity scale factor is obtained by a least- 
squares fit to the line of sight velocity dispersion of the 
model to Ab undances and Ra dial velocity Galactic Ori¬ 
gins Survey (|Ness et al.l 1201311 data for all stars within 
Galactocentric radius Rgc < 3.5 kpc at 6 = 5°, 7.5° and 
10° across |/| < 15°. We obtain a velocity scaling factor 
of 0.48. While these scalings lead to a model of roughly 
the right size and rotational velocity we stress that the 
model still does not match the MW and we only use it 
to qualitatively predict the expected trends in the MW, 
not their magnitude or precise location. 

We assume that the Sun is 8 kpc from the Galactic 
Genter, and place the observer at y = —8 kpc. We ori¬ 
ent t he bar at 27° to the line of sight (iWegg fc GerhardI 

Since we compare our model with APOGEE 
(jAlam et al.l 1201^ data, which targets bright red giant 
stars, we adopt a uniform selection function for star par¬ 
ticles at 2 kpc < Rit < 10 kpc, where R>, is the dis¬ 
tance from the Sun (jSchultheis et al.ir2014t iHavden et al.l 
l2015fl . Reducing the maximum R^ to 8 kpc does not sig¬ 
nificantly alter our conclusions. We use an opening angle 
of 0.5° for each LOSVD, to match the size of the smallest 
APOGEE bulge fields. The (off-plane) line of sight with 
the least particles contains over 2800 star particles while 
the best sampled (mid-plane) field has over 57,000 star 
particles; thus the shapes of the model LOSVDs are well 
determined. The top row of Eigure[T] shows the model’s 
surface density distribution. 

2.2. Line of Sight Velocity Distributions 

Viewing the model from the Solar perspective, we mea¬ 
sure the distribution of line of sight velocities in the 
Galactocentric restframe, Vgsr- Figure [T] shows the 
LOSVDs for various lines-of-sight (indicated in the top 
row) in the mid-plane (6 = 0 °, second row) and off-plane 
(b = 2°, third row). At ti each LOSVD at I < 12° has a 
single peak, both in the mid-plane and off the pl ane. The 
LOSV Ds have a shoulder to high Vgsr, which Hd et al.l 
(|2014ll showed is produced by stars at large distances 
seen close to tangentially. The peak in Vgsr moves to 
larger velocities with increasing /, but remains well be¬ 
low the Galaxy’s circular velocity. By 62 the LOSVDs at 


I = 8° and I = 10° have developed a second, high- Vgsr 
peak. This peak is more prominent than the Iow-Vgsr 
peak, due to the model’s very vigorous star formation in 
the nuclear disk, roughly ten times higher than in the 
MW for the corresponding region. This very high star 
formation rate quickly leads to a relatively massive nu¬ 
clear disk; thus the relative amplitudes of the low- and 
high-VcsR peaks are not predictions of the model. In¬ 
deed if we reduce the weight of star particles younger 
than 1 Gyr by a factor of 5, to compensate for the high 
star formation rate of the model, then the high-VcsR 
peaks become smaller than the main peaks, as seen in 
Figure [TJ The distribution around the high-VcsR peak 
is narrower {i.e. cooler) than that around the main peak 
and is skewed toward low Vgsr- Interior to 1 = 8°, the 
LOSVDs are broadened relative to those at ti, but no 
high- Vgsr peak is evident. At I > 14° no high- Vgsr peak 
is present in the mid-plane, indicating that the struc¬ 
ture responsible for the feature does not extend this far. 
The off-plane and mid-plane LOSVDs are not substan¬ 
tially different at ti, aside from the mid-plane hosting 
more stars at Vgsr > 100 kms“^. At 62 the high- Vgsr 
peaks, which dominate the mid-plane, are entirely ab¬ 
sent in the 6 = 2° LOSVDs. Therefore the presence of a 
nuclear disk is only evident in the mid-plane. As in the 
MW, outside the nuclear disk, the off-plane LOSVDs at 
{l,b) = (14°, 2°) also contain a statistically significant 
high- Vgsr peak/shoulder, but this is also present at ti, 
and is not related to the nuclear disk. Thus the kinematic 
signatures of a nuclear disk are (1) a second, high- Vgsr 
peak at roughly the circular velocity, (2) which is ab¬ 
sent a few degrees off the mid-plane, (3) is kinematically 
cooler than the Iow-Vgsr peak, and (4) is skewed toward 
low Vgsr- 

2.3. LOSVD Stacking 

The top row of Figure [1] shows color-coded maps of 
the average Vgsr, (Vgsr); the peak velocities at orbit 
tangent points manifest as the characteristic “winged” 
pattern of the (Vgsr) fields. Although the two (Vqsr) 
maps show the model before and after the nuclear disk 
forms, they are not very different, indicating that the 
formation of the nuclear disk does not lead to a wholesale 
change of the galaxy as much as populating new parts of 
its phase space. At the low longitudes of the nuclear disk, 
large (Vqsr) occurs only close to the galactic center while 
at other radii (Vgsr) is smaller. 

Even with a survey the size of APOGEE, the number 
of stars in individual fields is still relatively small, giv¬ 
ing a low signal-to-noise ratio for any second peak in any 
one field (|Li et al.l 1201411 . In order to overcome this dif¬ 
ficulty, we note that the Vgsr of the second peak does 
not change significantly with longitude at 6° < / < 10°. 
Therefore by stacking the LOSVDs we can enhance the 
signal-to-noise ratio of the high- Vgsr peak. Because the 
main peak is dominated by stars streaming along the 
bar, and (Vqsr) of these changes with I, the main peak 
in a stacked LOSVD will be quite broad. If we include 
Z < 4°, then the exponentially higher density of disk and 
bar stars near the center masks out any features at high 
Vgsr- In tfie bottom panels of Figure [T] we present a 
stack of the model’s LOSVDs at Z = 6°, 8° and 10°. 
As with the individual LOSVDs, a peak at high Vgsr is 
evident at O in the mid-plane but is absent at 6 = 2°. 





























y [kpc] 


Debattista et al. 


3 



Vgsr [km s '] [km s '] 

Fig. 1.— Top row: Face-on views of the model: contours indicate the surface density while colors show (Vgsr)- The bold dotted circles 
indicate the radii between which star particles are chosen (the selection function). The dashed lines show longitudes 4°-14° in 2° steps, 
color-coded as in the next two rows. Second row: Mid-plane LOSVDs for the different longitudes. Third row: LOSVDs at 6 = 2°, colored 
as in the rows above. Bottom row: Stacked model LOSVDs from Z = 6°, / = 8° and I = 10° in the mid-plane (black) and at & = 2° (red). 
The dashed black lines show the effect of reducing the weights of star particles younger than 1 Gyr by a factor of 5, to compensate for the 
very high star formation rate in the model. All LOSVDs have been normalized to unit peak. The left panels are at ti while the right ones 
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Field 

in 

bn 

N, 

Stack 

4336 

6.0 

0.0 

471 

mid-plane 

4355 

8.0 

0.0 

146 

mid-plane 

4365 

5.7 

2.0 

387 

off-plane 

4366 

5.7 

-2.0 

424 

off-plane 

4373 

7.8 

-2.0 

154 

off-plane 

4377 

7.7 

2.0 

149 

off-plane 


TABLE 1 

APOGEE FIELDS USED TO CONSTRUCT THE MID-PLANE AND 
OFF-PLANE STACKS. I AND b ARE THE GALACTIC LONGITUDE AND 
LATITUDE, RESPECTIVELY, OE THE FIELD CENTERS. N, IS THE 
NUMBER OF STARS SELECTED IN EACH FIELD. 



Eig. 2.— The stacked APOGEE LOSVDs for the mid-plane (red 
histogram) and at |fe| = 2° (blue histogram). The fields used are 
listed in Table 1. The black lines show the two Gaussians fitted to 
the mid-plane LOSVD. 

Moreover this second peak is still cooler than the low- 
Vgsr peak, and remains skewed toward it. Thus stacking 
LOSVDs preserves the kinematic signatures of a nuclear 
disk, and provides a reliable method for searching for a 
nuclear disk in the APOGEE data. 

3. APOGEE DATA 
3.1. Data Selection 

We select APOGEE survey stars in the fields of in¬ 
terest, excluding stars with the STAR BAD flag (corre¬ 
sponding to poor stellar parameter fits) and those flagged 
as flux and telluric standards. Stars with a velocity scat¬ 
ter between different visits of more than 5 kms“^ are also 
removed. (The same analysis including also stars flagged 
as STAR BAD, which leads to 763 in the plane and 1401 
out of the plane, gives results in agreement with the more 
conservative cut.) 

The small numbers of stars in the APOGEE commis¬ 
sioning data resulted in peaks with low signal-to-noise 
ratio. We increase the statistical significance of a high- 
bosR peak by stacking the APOGEE DR12 data in the 
longitude range 6° < I < 8° for fields in the mid-plane 
and off-plane at | 6 | = 2° (totalling 617 and 1114 stars, 
respectively). Table [T] lists the fields stacked together 
and the number of stars used from each held. 

Figure [5] plots these two stacked APOGEE LOSVDs. 
The mid-plane stack has a clear second peak at Vqsr 
220 kms“^, corresponding to rou ghly the circular v eloc- 
ity of the MW in the bulge region ( Sofue et al.ll2009ll . No 


comparable second peak is visible in the off-plane stacked 
LOSVD, which is non-Gaussian and ske wed toward hig h 
kbsR, *-e. it has a shoulder to high Vqsr (|Li et al.ll2014l) . 
A Kolmogorov-Smirnov test shows that the null hypoth¬ 
esis that the mid-plane and off-plane LOSVDs are drawn 
from the same distribution has a relatively low p-value 
of 0.04. 

We Ht two Gaussians to the mid-plane stacked LOSVD 
in the range —300kms“^ < Vgsr < 300kms“^, 
constrained such that the smaller Gaussian contains 
less than 25% of the stars (to avoid fitting just the 
skewed low-V qsr distribution with two Gaussians). We 
obtain a Iow-Vgsr component having mean velocity 
(Vgsr) = 24kms“^ and standard deviation ugsr = 
57 kms“^, while the high-V qsr component has (Vqsr) = 
217kms“^ and ugsr = 44kms“^, making it cooler 
than the Iow-Vgsr component. These two Gaussians 
are also shown in Figure [5] The velocity distribution at 
Vgsr > 200 kms“^ hints at a skewness opposite to that 
of the main distribution, but the signal-to-noise ratio is 
still too low for a robust measurement. 

The high-VcsR Gaussian has a significant number of 
stars associated with it, and is significantly separated 
from the Iow-Vgsr Gaussian. In order to test the like¬ 
lihood of such a second peak arising purely from Pois¬ 
son noise, we perform Monte-Carlo tests drawing 617 
stars from the off-plane stacked LOSVD. Fitting two 
Gaussians as before to the resulting LOSVD, we label 
as Gi and Gh the low- and high-VcsR components, re¬ 
spectively. We repeat this procedure 100,000 times, and 
for each we compute Nh/Ntot, the ratio of stars in the 
high-V gsr component to the total number of stars, and 
the overlap of the two components, defined as 

0 = J GiGhdVGSK- ( 1 ) 

The results are presented in Figure O the observed 
mid-plane stacked LOSVD has Nh/Ntot = 0.12 and 
O = 4.3. Only 0.025% of the Monte-Carlo samples 
have Nh/Ntot < 0 . 12 , while none of them have over¬ 
lap O < 8, showing that the observed double-peaked 
mid-plane stacked LOSVD is highly unlikely to result 
from Poisson noise. The APOGEE data therefore show 
a statistically significant double-Gaussian LOSVD in the 
mid-plane, the properties of which agree with 3 of the 4 
kinematic signatures of a nuclear disk from the simula¬ 
tion. While the signal-to-noise is too low to be sure if 
the high-V gsr peak is skewed to low Vgsr, the data are 
suggestive that it is. Therefore a kiloparsec-scale nuclear 
disk can explain the high-VcsR peaks in the APOGEE 
data. 

A simple estimate for the nuclear disk mass can be ob¬ 
tained from the fraction of stars in the high- Vgsr compo¬ 
nent of the double-Gaussian fit to the mid-plane LOSVD. 
If we conservatively assume that the nuclear disk mass 
contained within \z\ < 150 pc and 4° < |/| < 8 ° is 
12% of the total mass of the Besangon Galaxy model 
(|Robin et al.ll20ll) within this volume we obtain a lower 
limit to the mass of the nuclear disk ~ 5.8 x 10^ M 0 . 

4. DISCUSSION 

Attempts to explain the high-Vd sB peak directly via 
collisionless bar simulations fails (|Nidever et al.l l2012t 
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Fig. 3. — Frequency distribution of properties of double-Gaussian 
fits for the off-plane APOGEE stacked LOSVD sub-sampled Ns = 
10^ times to 628 stars. The side-panels indicate the distributions 
over the individual variables, normalized to unit peak. The param¬ 
eters for the fit to the mid-plane APOGEE stack are indicated by 
the filled red circle. In the side panels, the vertical dashed red lines 
indicate the values of Nh/Ntot and O for the mid-plane stack. 


iLi et al.l I2014D . However iMollov et al.l (|2015D demon¬ 
strated that resonant, bar-supporting 2:1 xl (with 
some mixture of higher order resonance) orbits by 
themselves can pro d uce s econd peaks. Subsequently 
lAumer fc Schdnrichl (|2015[ 1 argued that the selection 
function of APOGEE favors young stars recently trapped 
into resonant orbits. Their interpretation requires that 
the stars in the high-VoSR peaks are younger. The other 


main 2:1 resonant orbit family of bars, the x2 family, 
is orientated perpendicular to the bar. This family is 
enerally very poorly populated in the absence o f gas 


__ po pr __ 

Snarke fc SeilwoodlllQSTtIPfenniger fc Friedlilll991h . but 
when gas is present it is driyen inwards b y the bar and 
settles into x2 orbits (jBinney et al.llToMh . The gas can 
then form stars and produce nuclear rings and disks. 
We propose that the high-VcsR peak corresponds to 
a kiloparsec-scale disk composed of stars on x2 orbits. 
These orbits are stable and therefore our model does not 
require that the stars in the high-VcSR peak are young. 

Nuclear disks are know n in m any external galaxies 
(iScorza fc yan de n B oschl 119981: iZasoy fc Moiseeyl 

1999t Pizzella et al.l l2002t lEnisellem et al.l 12004 
Krainoyic et al.l 1200^ iLedo et al.l l20icilb the pres¬ 


ence of one in the MW is therefore not unusual. Nor is 
the kiloparsec scale unusual as a fraction of the bar size. 
For instance in NGC 3945 the ratio of se mi-major axes of 
the nuclear disk to bar is ^ 0.15 — 0.18 (lErwin fc Sparkd 
119991: iGole et 311120141) . whereas for the M W this ratio 
is ~ 0.2, if we adopt iWegg etldl (l2015D ’s 5 kpc bar. 
The gas ring in the simulation is ^ 5x larger than 
the MW’s CMZ, which i s coincident with a ste llar disk 
(|Launhardt et al.l 120021 : iSchonrich et al.l 120151 1. The 
large size of the gas ring in the model is a consequence 
of the still low resolution (50 pc) of our simulation 
(ILi et all i ormani et al.l 1201^ . This difference 

implies that the nuclear disk in the MW is not currently 
forming stars across its full extent. 

We anticipate that this proposal will inspire further de¬ 
tailed mapping of the central mid-plane of the MW. We 
will proyide predictions from our model of a kiloparsec- 
scale nuclear disk elsewhere. 
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